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ABSTRACT 

We at first introduce a novel photometric system, which consists of a Schmidt 
telescope, an objective prism, a CCD camera, and, especially, a multi-peak inter- 
ference filter. The multi-peak interference filter enables light in multi passbands 
to pass through it simultaneously. The light in different passbands is differen- 
tially refracted by the objective prism and is focused on the CCD separately, so 
we have multi "images" for each object on the CCD frames. This system enables 
us to monitor blazars exactly simultaneously in multi wavebands on a single tele- 
scope, and to accurately trace the color change during the variation. We used this 
novel system to monitor the BL Lacertae object S5 0716+714 during 2006 Jan- 
uary and February and achieved a very high temporal resolution. The object was 
very bright and very active during this period. Two strong flares were observed, 
with variation amplitudes of about 0.8 and 0.6 mags in the V band, respec- 
tively. Strong bluer- when-brighter correlations were found for both internight 
and intranight variations. No apparent time lag was observed between the V- 
and .R'-band variations, and the observed bluer-when-brighter chromatism may 
be mainly attributed to the larger variation amplitude at shorter wavelength. In 
addition to the bluer-when-brighter trend, the object also showed a bluer color 
when it was more active. The observed variability and its color behaviors are 
consistent with the shock-in-jet model. 

Subject headings: BL Lacertae objects: individual (S5 0716+714) — galaxies: 
active — galaxies: photometry 



1. INTRODUCTION 



Among the family of active galactic nuclei (AGNs), blazars manifest their prominent 
characteristic as rapid and strong variability nearly throughout the entire electromagnetic 
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wave bands. In the unified scheme of AGNs, blazars are those objec ts with their relativistic 



jets pointed basically to the observers (e.g., lUrry &: Padovanilll995l ). The jet is believed to 



originate from and be accelerated by a rotating supermassive black hole surrounded by an 
accretion disk. There are two types of blazars. Depending on whether or not they show 
strong emission lines in their spectra, they are called flat-spectrum radio quasars (FSRQs) 
or BL Lac objects. 

Several mechanisms, still in competition an d modification, haye bee n proposed to ex- 



Wagner fc Witzel 



Marscher fc Gear 



1 995]) . The most com 



19851 : bian et al.l fl~99lL 



plain the variability of blazars (see a review by 
mon suggestion is the shock-in-jet model (e.g 
in which shocks form at the base of the jet and propagate downstream, accelerating elec- 
trons and compressing magnetic fields, and resulting in the observed variabilit y. Alternative 
but less accepted scenar ios include the interstella r scintillation (ISS) (e.g.. iFiedler et al. 



19871 : iRickett et aDboOll ). the microle nsing effects dNottale!ll986l: ISchneider fc Weissl 



the instability in accretio n disk (e.g., lAbramowicz et al 



1991 



Chakrabarti &: Wiita 



98 



7) 



1993 



Mangalam fc Wiital Il993l ). and the geometrical effects (e.g., ICamenzind &: Krockenberger 



1992) . Some auth ors even suggested a combination among them (e.g., 



Wu et al. 



2005 



Dolcini et al.l 120051 ). For some special blazars, for instance, the BL Lac object OJ 287, 
which is th e only blazar that sho ws convincing evidences for a 12-year period in its optical 



variability (ISillanpaa et al 



1988) , binary black hole model was introduced to explain the 



observed periodicity (e.g., ISillanpaa et al.lll988l : IValtaoja et al.l l2000l ; iLiu fc Wul l2002l ). In 
brief, for the basic question in blazar studies, i.e., the variation mechanism, there exist vari- 
ous scenarios that still need to be verified and/or modified by future elaborate observations. 

In fact, an important but simple factor can discriminate among these various scenarios. 
It is the spectral in dex or color. The shock-in-jet mod el predicts spectral changes during 



the variations (e.g., iKirk. Rieger. fc Mastichiadid Il998l ). whereas the geometrical or light 



house effects do not. The ISS is frequency- dependent, but the micro-lensing effects are not. 
Moreover, these different mechanisms lead to different timescales of variability and different 
profiles of light curves. For example, the shock-in-jet and ISS will result in irregular light 
curves, the geometrical effects are usually linked to periodic variations, while micro-lensing 
effects should produce strictly symmetric light curves. By combining the spectral index, 
timescale, and light curve profile, the variation mechanism can be constrained. 

Therefore, in order to constrain the variability mechanisms in blazars, the crucial is to 
obtain the accurate spectral index. The spectral index is calculated by using two photomet- 
ric measurements in two different wavelengths. The traditional blazar monitoring on one 
telescope is to make observations in two or more wavelengths one by one, i.e., one can make 
an exposure with filter A at first, then change to filter B and make the second exposure. 
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After that one again change to filter C and make the third exposure, etc. The filters are 
usually used in a cyclic pattern. With this approach, one can achieve quasi-simultaneous 
observations in two or more wavelengths on a telescope. If the brightness of the target 
and the status of the instruments are both stable, one can reasonably calculate the spectral 
index with these quasi-simultaneous photometries. For blazars, however, the case is very 
different. Blazars are highly variable, with timescale as short as hours or even minutes. 
When the exposure with filter A is completed and a new exposure with filter B begins, the 
intrinsic brightness of the blazar may have changed. It will naturally result in large error 
or completely wrong result when calculating the spectral index by using the results of these 
two exposures. Therefore, all previous investigations related to spectral index or color of 
those highly active blazars, if based on quasi-simultaneous multi-waveband monitoring, have 
relatively low confidence level. 

In this paper, we introduce a novel photometric system, which includes a Schmidt 
telescope, an objective prism, a CCD camera, and a multi-peak interference filter. This 
new system enables us to make photometries exactly simultaneously in multi wavebands 
with a single exposure on a telescope, thus enables us to obtain the exact spectral index of 
blazars during the variations. The simultaneous exposures in multi wavebands also have the 
advantage of being invulnerable to the possible instability of instruments and weather. 

The BL Lac object S5 0716+714 is chosen as the first target of this new photometric 
system. It is very brig ht among all BL Lac objects, and it has a duty cycle of 1 (e.g., 



Wagner fc Witzell 119951). which m eans it is almost always active. In their monitorings of 



S5 0716+714, IVillata et al.l (120001 ) have determined a steepest recurrent variation slope as 



0.002 mag min 1 on both rising and decreasing phases, and iNesci. Massaro. fc Montagni 



(120021 ) found typical variation rate of 0.02 maghr 1 and a maximum r ising rate of 0.1 6 



maghr 1 . A maximum rising rate o f 0.1 maghr 1 was also reported by IWu et al.l (120051 ). 



Most recently, iMontagni et al.l (120061 ) claimed a highest rate of magnitude variation of 0.10 — 
0.12 mag h _1 . So this object with rapid and large-amplitude variation is very suitable to be 
monitored with our new photometric system. 

We monitored S5 0716+714 with this system in 2006 January and February. It is, to our 
knowledge, the first time that a multi-peak interference filter be used in blazar monitoring. 
Here we present the results, and based on accurate color tracing during the variation, we 
discuss the variation mechanism of this object. 
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2. ObSERVATIONS AND DATA REDUCTIONS 

2.1. The Photometric System 

The new photometric system consists of a Schmidt telescope, an objective prism, a CCD 
camera, and a multi-peak interference filter. The Schmidt telescope is 60/90 cm in diameter 
and is located at the Xinglong Station of the National Astronomical Observatories of China 
(NAOC). The CCD is a thick CCD, so it is not sensitive in blue wavelengths. It has a pixel 
size of 15 fim and a field of view of 58' x 58', resulting in a resolution of 1"7 pixel -1 . When 
this system is used for blazar monitoring, we read out only the central 512 x 512 pixels or 
14(5 x 14(5, which is adequate to cover the target blazar and the reference stars around it. 

The multi-peak interference filter is new in blazar monitoring and made by Changchun 
Institute of Optics, Fine Mechanics, and Physics, Chinese Academy of Sciences. The trans- 
mission curve of the filter is shown in Figure [TJ The left three passbands have narrower 
FWHMs than but similar central wavelengths to the traditional broad B, V, and R bands. 
Here we designate them as B', V, and R', respectively. There is still another passband 
beyond 9400 A. Because the CCD has a very low response at that wavelengths, data taken 
in this passband are of low quality and may only useful for bright sources. Table 1 gives the 
central wavelengths, maximum transmission, and bandwidths of the left three passbands. 

2.2. The Observations 

Figure [2] shows the finding chart of S5 0716+714 and an example CCD frame taken with 
this novel photometric system. Each source now has four "images" (actually four segments 
in the spectrum of the objective prism) on the frame. The four images, from bottom to 
top, correspond to the four peaks, from left to right, in Figure [H Because of the lower CCD 
responses and filter transmissions at the first and fourth peaks, the first and the fourth images 
are much fainter than the second and third ones. For the faint objects, only the second and 
third images are visible. Therefore, in the current paper, we do not use the data in the first 
and fourth passbands in the scientific analyses, but show two example light curves in the 
blue band. In the near future, a new thin CCD will be mounted to the telescope. It has 
much higher responses in the blue wavelengths, so we will have simultaneous photometries 
in at lease three passbands then. 

The monitoring of S5 0716+714 with this novel photometric system covered the period 
from 2006 January 1 to February 1. As a result of weather conditions, there are actually 
19 nights' data in total. The exposure time ranges from 120 to 540 s, depending on the 
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weather and moon phase. A histogram of the exposure time is shown in Figure [3j One can 
see that the majority of the exposure times do not exceed 300 s. The readout time is about 
5.6 s. So we have achieved a very high temporal resolution in multi wavebands with this 
new photometric system. The observational log and parameters are presented in Table 2. 



2.3. Data Reduction Procedures 

The procedures of data reduction include positional calibration, bias subtraction, flat- 
fielding, extraction of instrumental aperture magnitude, and flux calibration. The positional 
calibration was made for the i?'-band image, and the "coordinates" of the V- and S'-band 
images were known by their relative positions to their corresponding f?'-band images. The 
"coordinates" in the V and B' bands are, of course, not the coordinates in the common sense 
but are only labels for us to find the images during the photometry. The bias subtraction 
and flat-fielding were done as for the normal CCD frames. 

All images are elongated in the declinational direction, especially for the S'-band images. 
So apertures with special shapes were defined for them and are illustrated in Figure HI 
The traditional circular aperture is cut into two semi-circles along the RA direction, and 
a rectangle is inset into them. An elongated aperture is thus obtained. For the V- and 
i?'-band images, the widths of the rectangles (or the diameters of the semi-circles) are both 
4 pixels, but the heights are 5 and 4 pixels, respectively. For the S'-band image, the width 
and height of the rectangle are 3 and 14 pixels, respectively. All images have a sky annulus 
definition with inner and outer radii as 7 and 10 pixels, respectively. In order to eliminate 
the possible contamination of the apertures to their background, areas within 6 pixels in RA 
to the aperture center are excluded, and only the shaded areas are taken to be the actual 
sky background. 



Villata et al.l (119981 ) presented eight reference stars around S5 0716+714, as labeled in 
the left panel of Figure [21 We used the first four to calibrate S5 0716+714. Their broadband 
B, V, and R magnitudes were used to mimic their standard B', V, and R' ones during the 
calibrations. The results are presented in Table 2. The columns are observation date and 
time (UT), Julian Date, exposure time, and B', V, and R' magnitudes, their errors, and 
the differential magnitudes of reference stars. The differential magnitude of reference stars is 
defined as the magnitude difference (nightly mean set to 0) between star 3 and the average 
of the first four. It can be taken as a measure of the accuracy of the observations. 



There are, of course, some systematic differences between our B', V, and R' system 
and the broadband B, V, and R system. To quantitatively assess the differences, we carried 
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out photometric measurements for seven reference stars (1, 2, 3, 4, 6, 7, 8) by using the 
same aperture and background definitions and by using the same four reference stars for 
flux calibrations, as exactly did for S5 0716+714. Star 5 was excluded because its R! image 
is partly overlapped with the B' image of the BL Lac and its B' image is very close to the 
V and R! images of star 30. The photometric measurements of the seven stars were carried 
out for 20 best-observed frames taken on the photometric night of JD 2,453,742. Then the 
average B f , V, and R! magnitudes were calculated for the seven stars and compared with 
their B, V, and R magnitudes. The standard deviations of 7 (B — S')'s, 7 (V — V")'s, and 
7 (R — i?')'s were computed, and the results are 0.048, 0.028, and 0.036 mags, respectively. 
This indicates that there are only small systematic difference between our B', V, and R' 
system and the broadband B, V, R system. 



3. RESULTS 
3.1. Light Curves 

The V- and f?'-band light curves of the whole monitoring period are shown in Figure [5J 
The object was very active during this period. Two strong flares can be seen at JD 2,453,742 
and 2,453,757 with amplitudes of about 0.8 and 0.6 mags, respectively. Additional flares 
might also occur at JD 2,453,738, 2,453,745, and/or 2,453,749. The latter three are less 
certain because of the lack of observations on one or both sides of them. In addition to the 
strong internight variations, one can see clear intranight variations as well. Except for those 
less sampling nights (JD 2,453,745, 2,453,760, and 2,453,768), most intranight amplitudes 
are larger than 0.1 mags and can be as large as 0.3 mags. 

The whole light curve can be divided into two halves. Statistically, the first half has 
larger amplitudes in both internight and intranight variations than the second half has. So 
the object was more active in the first half than in the second. 

Figure [6] gives the intranight light curves on JD 2,453,737 and 2,453,742, which are 
respectively the first and the brightest nights. From top to bottom, they are in B', V, and 
R' bands, respectively. The large panels are the light curves of S5 0716+714, while the small 
ones give the differential magnitudes (nightly mean set to 0) between star 3 and the average 
of the first four. 



1 This is another source of error in the B' band data of S5 0716+714. Star 3 was not excluded just because 
its V and R' images are much brighter than the B 1 image of star 5 and thus are much less polluted. 
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On JD 2,453,737, the object kept on brightening, with some oscillations, then got fainter 
at the end. There was a sharp increase in brightness around JD 2,453,737.26. The V 
magnitude increased from 12.949 at JD 2,453,737.26147 to 12.886 at JD 2,453,737.27124 
(see Table 2), resulting in a brightening rate of about 0.004 magmin -1 , which is twice as 
fast as the rate reported by IVillata et al.l (120001 ) . On JD 2,453,742, the object underwent 
some small-amplitude oscillations at first. Then a flare was observed with an amplitude of 
about 0.235 mags in the V band and a duration of about 0.2 days. The flare peaks at 
JD 2,453,742.31079. It is also the peak of the first strong flare that lasted for five days (see 

Fig. ED. 

As can be seen, our temporal resolution is very high. The typical sampling intervals 
are 3.5 ~ 4.0 minutes in these two nights. On both nights, the light curves in the three 
bands are consistent excellently with each other, and the differential magnitudes in the V 
and R' bands have small rms's (0.0089 and 0.0093 for JD 2,453,737, and 0.0088 and 0.0103 
for JD 2,453,742), which both signify the accuracy of our observations. As expected, the 
S'-band data have large errors (for clarity, the error bars are not plotted), and the curves 
show frequent, irregular jumps in both large and small panels. So we did not show the B' 
band light curve in Figure El and will not use the B' band data in the following analyses. 



3.2. Correlations and Time Lags 



In order to investigate the correlation between the variations in different wave bands 
and to derive the possible lags between them, we at first carrie d out analyses of the z- 
transformed discrete correlation function (ZDCF; Alexander! 119971 ) on the V- and i?'-band 
variations. The ZDCFs between the V- and i?'-band variations are presented in Figure [7] for 
JD 2,453,737 and 2,453,742. The dashed lines are Gaussian fits to the points, and the dotted 
lines label the centers of the Gaussian profiles, which signify the time lags of the f?'-band to 
the l/'-band variations. The lags are 1.27 ± 0.22 and 0.06 ± 0.18 minutes for JD 2,453,737 
and 2,453,742, respectively. 



At the same time, the interpolated cross-correlation function (ICCF; lGaskell fc Peterson 



19871 ) was used to measure the time lags and the errors. The error was estimated with a 
model-independent Monte Carlo method, and the lag was taken as the centroid of the cross- 
correlation f unctions that were obta i ned with a large number of independent Monte Carlo 
realizations (IWhite &: Petersonlll994l ; iPeterson et al.lll998ll2004 ). Four thousand realizations 
were performed for both nights and gave the lags as 2.34±5.27 on JD 2,453,737 and — 0.01 ± 
1.88 minutes on JD 2,453,742. 
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All time lags are very short, shorter than our typical sampling interval, and are associ- 
ated with large errors. So we conclude that our observations do not detect apparent time 
lag between the V- and i?'-band variations. 

In optical regimes, variations at l ong wavelengths a r e usua lly reported to lag those at 
short wavelengths. For S5 0716+714, iQian. Tao. fc Fan! (120001 ) determined a n upper limit 



of lag of about 6 minutes between the V and / band variations. Similarly, IVillata et al. 



( 120001 ) claimed an upper limit of 10 minutes to the possible delay betw een the B and I band 
variations using high quality data. Most recently, IStalin et al.l (120061 ) reported time lags of 
about 6 and 13 minutes for the V and R band variations on two individual nights. But they 
also noted that their measurement intervals were close to these putative lags, and the lags 
should be treated with caution. For other objects, for example, iRomero. Cellone. Combi 
( I2000l ) detected time lags of a few to 17.3 minutes between the V and R band variations 
for PKS 0537—441. However, the monitoring on 3C 66A does not strongly support any 
time lag between variations at different wavelengths. The detailed spectral modeling with a 
time-dependent leptonic one-zone jet model also does not predict such spectral hysteresis (M. 
Bottcher, private communication). Our data, characterized by very high temporal resolution 
and simultaneous measurements in multi frequencies, do not detect apparent lag between 
variations in the V and R' bands. Compared to the much longer lags between variations in 
X-rays (e.g.. iTakahashi et al.lll996l . ~ 1 hr) or between X-ray and EUV/UV variations (e.g., 



Urry et al.l 1 199 71 . 1 ~ 2 days), these very short or even absent time lags in optical regimes 
may be the result of very small frequency intervals, and may indicate that the photons in 
these wavelengths should be produced by the same physical process and emitted from the 
same spatial region. 



3.3. Color Behavior 

3.3.1. The Color- Magnitude Correlation 

As mentioned in §1, the color behavior during the variations can put strong constraints 
on the variation mechanisms. Figure [S] displays the color versus magnitude distributions 
on JD 2,453,737 and 2,453,742. There is a clear bluer-when-brighter (BWB) chromatism 
on both nights. The dashed lines are linear fits to the points. The Pearson correlation 
coefficients are 0.368 and 0.444, and the significance levels are 1.1 x 10~ 6 and 1.6 x 10~ 8 , 
respectively, which indicate strong correlations between color and magnitude. The dashed 
lines have slopes of 0.060 and 0.076 for JD 2,453,737 and 2,453,742, respectively. 

For all nights, there is also a BWB trend for S5 0716+714. This is illustrated in 
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Figure [9l The dashed line is the linear fit to the 1815 points, and it has a slope of 0.061. 
The Pearson correlation coefficient is 0.552, indicating a strong correlation between the color 
and magnitude. 

The color or spectral behaviors of S5 0716+714 and other blazars have been investi- 



gated in optica l bands by many authors (IVagnetti. Trevese. fc Nescill2003l ; IStalin et al.ll2006 



Wu et al.ll2005l . and references therein). On short timescales, the blazars usually BWB chro- 
matism when they are in an active or flaring state. On long timescales, an achromatic trend 
is usually found. Our intranight color behavior is in agreement with previous results, whereas 
the long-term internight variations show a different color behavior. Our BWB chromatisms 
argue strongly against the microlensing or geometrical effects as the variation mechanisms 
in S5 0716+714. 

In the variability of blazars, two factors may result in color change. One is the difference 
in variation steps at different frequencies, and the other, the difference in variation amplitudes 
of these variations. If the variation at high frequency leads that at low frequency, we will 
observe a BWB chromatism. On the other hand, if the fluxes at different frequencies vary 
simultaneously, but the amplitude is larger at high frequency than at low one, we will again 
observe a BWB chromatism. In the actual cases, it may be a combination of these two 
factors that results in the color changes of blazars. In §3.2, we do not find apparent time 
lag between variations in different wave bands. The observed BWB phenomena should be 
mainly attributed to the larger amplitudes at higher frequency This will be confirmed in the 
next section. 



3.3.2. Color- Activity Correlation 

The bottom panel of Figure [5] displays the temporal evolution of the color V — R! or 
the V — B! "color curve". Most nights have a scatter less than 0.1 in color, except for 
JD 2,453,749 and 2,453,754. The two nights are close to or within the phase of full moon, so 
the photometries have relatively larger errors, and the colors are much more scattered than 
those in other nights. 

From the color curve, one can see that the object generally showed a bluer color when 
it was brighter, which is consistent with the result in Figure IHJ However, it is also clear that 
the color does not correlate solely with the brightness. For example, the object was brightest 
but not bluest on JD 2,453,742. The brightness on JD 2,453,737, 2,453,738, and 2,453,740 
is much fainter than that on JD 2,453,742, but the former nights have comparable or even 
bluer colors than the latter night has. In analogy to the division of the light curve, the color 
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curve can also be divided into two halves, i.e., those before and after JD 2,453,750. The first 
half is much bluer than the second, as manifested by the dashed line, the mean color, in the 
color curve in the bottom panel of Figure [5J Therefore, the object seems show a bluer color 
when it is more active. 

To study the relation between the color and the activity in more detail, we at first 
need to define the activity of a blazar, or how a blazar is called active. The activity has 
two aspects; one is the frequency or rate, and the other, the amplitude. The frequency of 
activity can be described by duty cycle, which is d efined as the fraction o f time a source 



spends more than 3a away from its weekly average (jWagner fc Witzelll!995l ) 



For the amplitude of activity, we define it as the sum of the intranight and internight 
amplitudes 

^4 ^intran ^intern- (1) 

The intranight amplitude can be simply calculated as the difference between the maximum 
and minimum magnitudes (or fluxes) of a certain nig ht, as did in §3.1. 3 . A m ore precise way 



is to take the measurement error into account. As in lHeidt fc Wagner! (119961 ). the intranight 
amplitude is defined as 

^■intran ^/(?77, max H-min)^ 2(7^, (2) 

where m max and m m \ n are respectively the maximum and minimum magnitudes (or fluxes) 
of a certain night, and a is the measurement error. Here we take a as the nightly average of 
the differential magnitude of reference stars (SV or 5R' in Table 2). 

For the internight amplitude, the average magnitudes (or fluxes) are calculated at first 
for all nights. Then the internight amplitude of the zth night is defined as 

A _ \( m )i ~ Mt-il + l( m )» ~ ( o\ 

-^■intern, i 2 ' \ ) 

where (ra)i_i, (m)i, and (m) i+ i are the average magnitudes of the (i — l)th, ith, and (i + l)th 
nights. When there is no observation on the (i — 1) or (i + 1) night, the data on the night 
closest to the (i — 1) or (i + 1) night are used if the time differences are not large, but the 
resulting amplitude may have large error. For the first and the last night of the sequence, 
only the first or second absolute value is calculated in Equation 3 and it is not divided by 2. 
In this paper, we will focus on the amplitude rather than the frequency of activity. So we 
call a blazar "active" when it has a large amplitude of activity. 

With these definitions, we calculated the intranight and internight amplitudes of ac- 
tivity of S5 0716+714 on the 19 monitoring nights and listed the results in Table 3. The 
columns are Julian Date, average V magnitude, internight and intranight amplitudes in the 
V band, average B! magnitude, internight and intranight amplitudes in the R' band, and 
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time duration of the monitoring in that night. Figure [10] illustrates the relation between the 
color and the amplitude of activity. We do not plot the data on JD 2,453,749 and 2,453,754 
because they have much larger measurement errors and much larger scatters in color than 
those on other nights. Data on JD 2,453,745, 2,453,760, and 2,453,768 are also excluded 
from the figure because they are much less sampled than those on other nights. So there are 
14 nights' data left. The triangles and circles designate nights of the first and second halves 
of the light curve, respectively. The former has much bluer color than the latter, as expected. 
The solid line is the best fit to the 14 points. The Pearson correlation coefficient is —0.851 
and the significance level is 1.15 x 10~ 4 , which mean strong correlation between color and 
amplitude of activity. The correlation confirms the result of the above visual inspection: the 
object showed a bluer color when it was more active. 

Figure [TT] displays whether the intranight and internight amplitudes separately correlate 
with the color. The distribution on the left panel has a Pearson correlation coefficient of 
—0.745 and a significance level of 2.24 x 10~ 3 , while the right panel shows a correlation 
coefficient of —0.576 and a significance level of 0.03. The former correlation is stronger than 
the latter, but they are both much weaker than the correlation in Figure [Ull 

Most observations detected a BWB chromatism in blazars. At the same t ime, a minor- 



ity o f monitorings fo und a redder- when-brighter trend in some blazars (e.g., iRaiteri et al. 



20031 ; iGu et al.ll2006l ). This inconsistency may be explained with the bluer- when-more-active 
(BWMA) chromatism: even though the object is in a faint state, it can still have relatively 
large amplitude of variations, so it can show a relatively blue color. 



4. COLOR EVOLUTION DURING THE FLARE 

According to theoretical predictions, the color of a blazar usually evolves with a wave- 
like pattern during a flare, and there is usually a loop path on the color versus magnitude 
(or flux) diagram, in either clockwise or anti-clockwise direction. The direction depends 
on the frequency that the observation is made and the peak frequency of the synchrotron 



Kirk. Rieeer. & Mastichiadis 1998). The loop path has been frequently observed in X-ray 


(e.g., Sembav et a 


1993: Takahashi et al.lll996: 


Kataoka et al. 2000: 


Zhansr et al. 1999. 2002: 


Malizia et alJl2000: 


Ravasio et all2004f). and. in one case, in infrared ( 


Gear, Robson, & Brown 


1986f). In optical regimes, only Xilouris et al. ( 


2006) have reported a similar pattern. 



In our monitoring, we recorded two strong flares, i.e., the one from JD 2,453,740 to 
2,453,744 and the other from JD 2,453,754 to 2,453,760 (see Fig. ED . The temporal evolutions 
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of the color do follow a wave-like pattern during the two flares (see the bottom panel of Fig. [5]). 
On the color versus magnitude diagrams in Figure [121 the results are controversial. When 
the individual measurements are considered, both diagrams are just messy; when we focus 
only on the nightly means, there is a loop in clockwise direction on both diagrams, as labeled 
by the arrows. 

The loop path can be explained by a spectral hysteresis, i.e., the variation in one band 
lags that in the other. Our observations do not find apparent time lags between the intranight 
variations in different wave bands, so the individual measurements do not show a loop path. 
However, the loop paths described by the nightly means may suggest that S5 0716+714 has a 
spectral hysteresis in its long-term (from hou rs to days) variations . By using t he same Monte 
Carlo algorithm of IWhite &: Peterson! (119941 ) and iPeterson et al.l (119981 . 12004 ) , we derived a 
time lag of 0.10 ± 3.98 days between the V- and i?'-band nightly-mean light curves. Given 
the large error and the short lag relative to the sampling interval (1 day), the result is of 
relatively low significance level in quantity, but should be qualitatively reasonable. 



5. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have introduced a novel photometric system, which consists of a 
Schmidt telescope, an objective prism, a CCD camera, and, especially a multi-peak inter- 
ference filter. The multi-peak interference filter enables light in multi passbands to pass 
through it simultaneously. The light in different passbands are differentially refracted by the 
objective prism and focus on the CCD separately, so we have multi images for each object 
on the CCD frames. This system enables us to carry out simultaneous photometric measure- 
ments in multi wavebands, and it has the advantage of being invulnerable to the possible 
instability of instruments and weather. When using for blazar monitorings, the system can 
accurately trace the color change during the variation, which is crucial in constraining the 
variation mechanisms of blazars. 

The multi-peak interference filter is new, at least in blazar monitorings. The transition 
from the traditional one-passband filter to the new multi-peak interference filter is analogous 
to the transition from the traditional slit spectrograph to the fiber spectrograph. The new 
filter can greatly increase the efficiency of data sampling and is worthy being extended to 
the observations of other targets, especially those variable sources. 

This novel photometric system was used in the monitoring of S5 0716+714 and we 
achieved a very high temporal resolution. The BL Lac object was very bright and very 
active during our monitoring period and showed strong variations on both internight and 
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intranight timescales. Two strong flares were observed, with amplitudes of variations of 
about 0.8 and 0.6 mags, respectively. There may be three more but less certain flares. 
Significant BWB correlations were found for both internight and intranight variations, and 
the BWB internight variability is diff erent from the basically achromatic long-t erm variability 
of S5 0716+714 observed previously (jGhisellini et al.lll997l ; iRaiteri et al.ll2003l ). No apparent 
time lag was found between the V- and i?'-band variations. The BWB chromatism should 
be mainly attributed to the larger amplitude of variation at higher frequency. 

In addition to the BWB chromatism, we also found a BWMA chromatism, which is new 
and has not been reported before. It helps to explain some exceptional color behaviors (e.g., 
redder-when-brighter) observed in blazars sometimes. The theoretical reason for the BWMA 
chromatism is still unknown. Phenomenologically, when the object is more active (has larger 
variation amplitude), the difference in variation amplitudes in different wavebands is larger, 
so the object is bluer. Unlike the BWB chromatism, which can be either long- or short-term 
behavior, the BWMA phenomenon is only a long- or, at least, intermediate-term behavior. 
To reveal or confirm the BWMA chromatism, the monitoring should cover a relatively long 
period and consists of as many successive nights as possible. Intensive intranight observations 
are also desired. 

The color behaviors can put strong constraints to the variation mecha nisms of blazars. 
The m echanisms can be broadly classified into extrinsic and intrinsic origins (IWagner &: Witzel 
19951 ). The extrinsic origins include the ISS and gravitational microlensing effects. The ISS 
can explain variations at low radio frequencies, but can not produce variations at optical 
wavelengths. The microlensing is an achromatic process and should result in strictly sym- 
metric light curves, which are both contrary to our results. So the variations of S5 0716+714 
observed by us are unlikely to have an extrinsic origin. The intrinsic origins include the 
shock-in-jet scenario, accretion disk instability, and geometrical effects. In the accretion disk 
scenario, external gravitational perturbations may induce flares and hot spots in the accre- 
tion disk a round a black hole. This may account for the optical-U V microvariability observed 
in blazars (IChakrabarti fc Wiitalll993l ; iMangalam &: Wiitalll993l ). but it has difficulty in ex- 
plaining the strong or rapid flares in blazars, as observed by us. I t also can not account fo r 
the correlated radio-optical variations observed in S5 0716+714 (IQuirrenbach et al.lll99ll ). 
In fact, the general absence of the big blue bump in the SED of blazars also argue against a 
significant contribution of the thermal emission from the accretion disk to the overall emis- 
sion of blazars. The geometrical effects are that the helical or precessing jet leads to varying 
Doppler boosting towards the observers and hence the flux variations. They usually result 
in achromatic and periodic variations, which is also not the case in our results. Now the 
other mechanisms are all excluded, the shock-in-jet scenario becomes most plausible. In fact, 
the irregular variations and, especially, the BWB color behaviors are all consistent with the 
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predictions of the shock-in-jet model. 

Variation mechanism is an essential issue in the investigations of blazars. Our new 
photometric system has the advantages of accurate color trace, high temporal resolution, 
and being invulnerable to instrumental and weather instabilities. A new thin CCD camera 
will soon be mounted to the telescope, and it has much higher responses in blue wavelengths 
than the current thick CCD has. The objective prism may also soon be replaced by an 
diffraction grating. The upgraded system will be used to monitor a few bright blazars in 
the near future, and the simultaneous three-band photometries will help to investigate the 
variation mechanisms in them with high confidence level. 

The authors thank the anomymous referee for constructive suggestions and insightful 
comments. We thank Dr. Bradley M. Peterson for sending the code for calculating the time 
lags between variations in different wave bands. This work has been supported by the Chinese 
National Natural Science Foundation grants 10603006, 10473012, 10573020, 10633020 and 
10303003. 
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Fig. 1. — Transmission curve of the multi-peak interference filter. 



Fig. 2. — Finding chart of S5 0716+714 (left) and example CCD frame (right) taken with 
the novel photometric system. The sizes are 10' x 10' and 14'5 x 14f5, respectively. Labeled 
on the finding chart are the BL Lac object and eight reference stars. The first four stars are 
used to calibrate the BL Lac. Each object has four images on the frame, as indicated for 
S5 0716+714. For those faint sources, only two central images are visible. 
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Fig. 3. — Histogram of exposure time. Most exposure times are equal to or less than 300 s. 
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Fig. 4. — Aperture and sky annulus definitions for an example image on a frame. The tradi- 
tional circular aperture is cut into two semi-circles along the RA direction, and a rectangle 
is inset into them. An elongated aperture is thus obtained. Images in different colors have 
different sizes of rectangles and semi-circles. The shaded areas in the annulus are taken to 
be the sky background. See text for details. 
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Fig. 5. — Light curves in the V and R! bands and color evolution. The first halves of the 
light curves have much larger amplitudes in both internight and intranight variations than 
the second halves. The colors in the first half is much bluer than those in the second, as 
manifested by the dashd line, the mean color. 
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Fig. 6.— Intranight light curves on JD 2,453,737 and 2,453,742, in B', V, and R' bands 
from top to bottom. The large panels show the light curves of S5 0716+714, while the small 
ones display the differential magnitudes between star 3 and the average of the first four. The 
errors are not plotted for clarity. The B' band light curves have large errors, as can be seen 
in the small panels. 
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Fig. 7.— ZDCFs between the V- and fl'-band variations on JD 2,453,737 and 2,453,742. 
The dashed lines are Gaussian fits to the points, and the dotted lines (at 1.27 and 0.06 
minutes for left and right, respectively) label the centers of the Gaussian profiles, which 
signify the lags of the R' band to V band variations. 




Fig. 8. — Color versus magnitude distribution for JD 2,453,737 and 2,453,742. The dashed 
lines are linear fits to the points. The Pearson correlation coefficients of 0.368 and 0.444 
indicate strong correlations between color and magnitude. The slopes of the two fits are 
0.060 and 0.076, respectively. 
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Fig. 9. — Color versus magnitude distribution for all nights. The linear fit to the points, 
the dashed line, gives a slope of 0.061. The Pearson correlation coefficient of 0.552 means a 
strong correlation between color and magnitude. 



- 26 - 





I I I I I | I I I 

r=-0.851 






0.42 




A , 


- 




- 


A a 


- 


0.44 


A 










A 






O 






0.46 


o ° 








- o ° 






0.48 


1 1 1 1 1 1 1 1 1 1 




1 



0.2 0.3 0.4 0.5 0.6 

A 



Fig. 10. — Correlation between color and amplitude of activity. The triangles are nights 
of the first half of our monitoring period, whereas the open circles are nights of the second 
half. The solid line is the best fit to all 14 points. The Pearson correlation coefficients, 
r = —0.851, indicate a strong trend that the object become bluer when it is more active. 
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Fig. 11. — Correlations between color and intranight (left) and internight (right) amplitudes. 
The former correlation is stronger than the latter, but both are much weaker than the 
correlation in Fig. [TU1 
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Fig. 12. — Color versus magnitude distributions during two flares. Different symbols signify 
measurements on different nights. The individual measurements do not describe a loop path, 
while the nightly means define two loop paths in clockwise direction, as indicated by the 
arrows. 
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Table I. Parameters of the Left Three Passbands of the Filter 



Passband Wavelength Max Trans FWHM 

(A) (%) (A) 



B' 


4575.2 


59.40 


355 


V 


5672.9 


88.28 


210 


B! 


6786.6 


58.34 


162 



Table 2. Observational log and results in the B', V, and R! bands. 



Date 


Time 


Julian Date 


Exp. 


B' a 


R' 

err 


SB 1 


V 


V 


SV 


Rf 


TV 

-TL err 


SR' 








(s) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


2006 01 01 


11:30:07 


2453736.97925 


240 


13.463 


0.027 


0.022 


13.081 


0.009 


0.009 


12.662 


0.009 


-0.014 


2006 01 01 


11:49:27 


2453736.99268 


180 


13.476 


0.027 


0.039 


13.072 


0.011 


0.019 


12.644 


0.012 


0.000 


2006 01 01 


11:54:06 


2453736.99585 


180 


13.521 


0.031 


0.061 


13.059 


0.010 


-0.018 


12.628 


0.011 


-0.006 


2006 01 01 


12:30:03 


2453737.02075 


180 


13.497 


0.026 


0.053 


13.051 


0.010 


-0.010 


12.597 


0.010 


-0.006 


2006 01 01 


12:33:27 


2453737.02319 


180 


13.531 


0.028 


-0.015 


13.045 


0.010 


0.000 


12.619 


0.009 


0.011 



Note. — Table 2 is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. 

a The B' band data have large errors and should be used with great caution. 
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Table 3. 


Average 


Magnitudes and Variation Amplitudes on 


Individual Nights 


Julian Date 


(V) 


-"■V, intern 


A a 
-^V'jintran 


(R>) 


^4r', intern 


A a 
^R^intran 


Duration 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(day) 


2453737 


12.952 


0.182 


0.302 


12.533 


0.186 


0.298 


0.444 


2453738 


12.770 


0.275 


0.199 


12.347 


0.263 


0.202 


0.403 


2453740 


13.138 


0.380 


0.098 


12.688 


0.350 


0.091 


0.317 


2453741 


12.746 


0.302 


0.269 


12.329 


0.289 


0.246 


0.423 


2453742 


12.533 


0.207 


0.278 


12.109 


0.201 


0.251 


0.421 


2453743 


12.733 


0.181 


0.095 


12.292 


0.168 


0.084 


0.411 


2453744 


12.895 


0.162 


0.174 


12.446 


0.154 


0.143 


0.429 


2453745 


12.691 


0.204 


0.055 


12.257 


0.189 


0.049 


0.076 


2453749 


13.062 


0.186 


0.242 


12.606 


0.179 


0.212 


0.223 


2453754 


13.063 


0.044 


0.077 


12.614 


0.053 


0.075 


0.359 


2453756 


12.977 


0.160 


0.151 


12.516 


0.150 


0.131 


0.196 


2453757 


12.817 


0.117 


0.191 


12.366 


0.110 


0.184 


0.402 


2453758 


12.892 


0.135 


0.103 


12.436 


0.125 


0.096 


0.420 


2453759 


13.087 


0.162 


0.142 


12.617 


0.160 


0.122 


0.334 


2453760 


13.171 


0.065 


0.055 


12.712 


0.070 


0.055 


0.087 


2453761 


13.216 


0.069 


0.179 


12.756 


0.074 


0.161 


0.417 


2453762 


13.226 


0.052 


0.160 


12.764 


0.047 


0.151 


0.405 


2453767 


13.320 


0.094 


0.116 


12.849 


0.085 


0.128 


0.409 


2453768 


13.387 


0.067 


0.065 


12.915 


0.066 


0.046 


0.026 



a The spurious measurements were excluded when calculating the intranight amplitudes. 



